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ABSTRACT 

We repeat the directional spherical real Morlet wavelet analysis, used to detect 
non-Gaussianity in the Wilkinson Microwave Anisotropy Probe (WMAP) 1-year data 
(McEwen et al. 2005a), on the WMAP 3-year data. The non-Gaussian signal previously 
detected is indeed present in the 3-year data, although the significance of the detection is 
reduced. Using our most conservative method for constructing significance measures, we find 
the significance of the detection of non-Gaussianity drops from 98.3 +0.4% to 94.9+0.7%; the 
significance drops from 99.3 + 0.3% to 97.2 + 0.5% using a method based on the;if^ statistic. 
The wavelet analysis allows us to localise most likely sources of non-Gaussianity on the sky. 
We detect very similar localised regions in the WMAP 1-year and 3-year data, although the 
regions extracted appear more pronounced in the 3-year data. When all localised regions are 
excluded from the analysis the 3-year data is consistent with Gaussianity. 
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1 INTRODUCTION 

Recent measurements of the cosmic microwave background 
(CMB) anisotropics, in particular those made by the Wilkinson 
Microwave Anisotropy Probe (WMAP), provide data of unprece- 
dented precision with which to study the origin of the universe. 
Such observations have lent strong support to the standard cos- 
mological concordance model. Nevertheless, many details and as- 
sumptions of the concordance model are still under close scrutiny. 
One of the most important and topical assumptions of the standard 
model is that of the statistics of the primordial fluctuations that give 
rise to the anisotropics of the CMB. In the simplest inflationary 
models, primordial perturbations seed Gaussian temperature fluc- 
tuations in the CMB that are statistically isotropic over the sky. 
However, this is not necessarily the case for non-standard inflation- 
ary models or various cosmic defect scenarios. 

The assumptions of Gaussianity and isotropy have been ques- 
tioned recently with many works highlighting deviations from 
Gaussianity in the WMAP 1-year data (WMAPl; Bennett et al. 
2003), calculating measures such as the bispectrum and Minkowski 
functionals (Komatsu et al. 2003; Magueijo & Medeiros 2004; 
Land & Magueijo 2005a; Medeiros & Contaldi 2005), the genus 
(Colley & Gott 2003; Eriksen et al. 2004), correlation functions 
(Gaztanaga & Wagg 2003; Eriksen et al. 2005; Tojeiro et al. 2006), 
low-multipole alignment statistics (de Oliveira-Costa et al. 2004; 
Copietal. 2004, 2005; Schwarz et al. 2004; Slosar & Seljak 
2004; Weeks 2004; Land & Magueijo 2005b,c,d,e; Bielewicz et al. 
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2005; de Oliveira-Costa & Tegmark 2006), structure alignment 
statistics (Wiaux et al. 2006) phase associations (Cliiang et al. 
2003; Chiang & Naselsky 2004; Coles et al. 2004; Dineen et al. 
2005), local curvature (Hansen et al. 2004; Cabella et al. 2005), 
the higher criticism statistic (Cayon et al. 2005), hot and cold 
spot statistics (Larson & Wandelt 2004, 2005), fractal statis- 
tics (Sadegh Movahed et al. 2006) and wavelet coefficient statis- 
tics (Vielva et al. 2003; Mukherjee & Wang 2004; McEwen et al. 
2005a, 2006; Cruz et al. 2005, 2006a). Some statistics show consis- 
tency with Gaussianity, whereas others provide some evidence for 
a non-Gaussian signal and/or an asymmetry between the northern 
and southern Galactic hemispheres. Although the recently released 
WMAP 3-year data (WMAP3; Hinshaw et al. 2006) is consistent 
with the WMAPl data, a more thorough treatment of beams, fore- 
grounds and systematics in the 3-year data, in addition to a fur- 
ther two years of observing time, mean that WMAP3 provides a 
more reliable data-set on which to confirm or refute previous re- 
sults. A Gaussianity analysis is performed on the WMAP3 data 
by Spcrgcl et al. (2006), using the one point distribution function, 
Minkowski functionals, the bispectrum and the trispectrum. No 
evidence is found for non-Gaussianity; however, the authors do 
not re-evaluate the large number of statistical tests that have been 
used to detect non-Gaussianity in the WMAPl data. Indeed, devi- 
ations from Gaussianity and isotropy have recently been detected 
in the WMAP3 data, using measures such as anisotropy statis- 
tics (Helling et al. 2006; Bernui et al. 2006) phase associations 
(Chiang et al. 2006) and wavelet coefficient statistics (Cruz ct al. 
2006b), with little change in the significance levels obtained by 
each technique for the first and third year data. Although the depar- 
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tures from Gaussianity and isotropy detected in the WMAPl and 
WMAP3 data may simply highlight unremoved foreground con- 
tamination or other systematics, which itself is of importance for 
cosmological inferences drawn from the data, if the source of these 
detections is of cosmological origin then this would have important 
implications for the standard cosmological model. 

In this letter we focus on the significant detection of non- 
Gaussianity that we made previously in the WMAPl data using 
directional spherical wavelets (McEwcn ct al. 2005a), to see if the 
detection is still present in the WMAP3 data. The remainder of 
this letter is organised as follows. In Section 2 we briefly review 
the analysis procedure and discuss the data maps considered. Re- 
sults are presented and discussed in Section 3, before concluding 
remarks are made in Section 4. 



2 NON-GAUSSIANITY ANALYSIS 

We repeat on the WMAP3 data our non-Gaussianity analysis per- 
formed previously on the WMAPl data (McEwcn ct al. 2005a), fo- 
cusing only on the most significant detection of non-Gaussianity 
made previously. We refer the reader to our previous work 
(McEwen et al. 2005a) for a detailed description of the analysis 
procedure and present here only a very brief overview. 

We apply a spherical wavelet analysis to probe the WMAP 
data for non-Gaussianity. Wavelets are an ideal tool to search for 
deviations from Gaussianity due to the scale and spatial localisa- 
tion inherent in a wavelet analysis. To perform a wavelet analy- 
sis of full-sky CMB maps we apply our fast continuous spherical 
wavelet transform (CSWT; McEwen et al. 2005b), which is based 
on the spherical wavelet transform developed by Antoine, Van- 
dergheynst and colleagues (Antoine & Vandergheynst 1998, 1999; 
Antoine et al. 2002, 2004; Wiaux et al. 2005) and the fast spher- 
ical convolution developed by Wandelt & Gorski (2001). We use 
only the real Morlet wavelet in this analysis since it gave the 
most significant detection of non-Gaussianity in the WMAPl data 
(McEwen et al. 2005a). 

To minimise the contribution of foregrounds and systemat- 
ics to CMB anisotropy measurements, the WMAP assembly con- 
tains a number of receivers that observe at a range of frequencies 
(Bennett et al. 2003). In this analysis we consider the signal-to- 
noise ratio enhanced co-added map constructed from the WMAP3 
data. This map is constructed by the same procedure described gen- 
erally by Komatsu ct al. (2003) and described in the context of our 
non-Gaussian analysis by McEwen et al. (2005a). We use the fore- 
ground reduced sky maps and apply the KpO mask to remove resid- 
ual Galactic emission and known point sources. The foreground 
maps and mask are available from the Legacy Archive for Mi- 
crowave Background Data Analysis (LAMBDA) website'. 

To quantify the significance of any deviations from Gaussian- 
ity we perform 1000 Monte Carlo simulations. This involves simu- 
lating 1000 Gaussian co-added maps. Each simulated map is con- 
structed in an analogous manner to the co-added map constructed 
from the data. A Gaussian CMB realisation is simulated from the 
theoretical power spectrum fitted by the WMAP team (the power 
spectrum we use is also available from LAMBDA). Measurements 
made by the various receivers are then simulated by convolving 
with realistic beams and adding anisotropic WMAP3 noise for 
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Figure 1. Real Moriet wavelet coefficient skewness statistics (y = 72°). 
Points are plotted for the WMAPl data (solid, green, squares), WMAP3 
data (solid, blue, circles) and the WMAP3 data with localised regions re- 
moved (dashed, blue, triangles). Confidence regions obtained from 1000 
WMAP3 Monte Carlo simulations are shown for 68% (red), 95% (orange) 
and 99% (yellow) levels, as is the mean (solid white line). 



each receiver. The simulated observations for each receiver are then 
combined to give a co-added map. 

To probe the WMAP3 data for deviations from Gaussianity the 
skewness of the real Morlet wavelet coefficients is examined over 
a range of scales and orientations (the scales and orientations con- 
sidered are defined in McEwen et al. 2005a). Any deviation from 
zero is an indication of non-Gaussianity in the data. An identical 
analysis is performed on the 1000 Gaussian simulations to quantify 
the significance of any deviations. 



3 RESULTS 

The skewness of the real Morlet wavelet coefficients of the co- 
added WMAP3 map are displayed in Figure 1, with confidence 
intervals constructed from the 1000 WMAP3 Monte Carlo simu- 
lations also shown. Only the plot corresponding to the orientation 
of the maximum deviation from Gaussianity is shown. The non- 
Gaussian signal present in the WMAPl data is clearly present in the 
WMAP3 data. In particular, the large deviation on scale an = 550' 
and orientation y = 72° is almost identical (although it is in fact 
very marginally lower in the WMAP3 data). 

Next we consider in more detail the most significant deviation 
from Gaussianity on scale an = 550' and orientation y = 72°. Fig- 
ure 2 shows histograms of this particular statistic constructed from 
the WMAPl and WMAP3 Monte Carlo simulations. The measured 
statistic for the WMAPl and WMAP3 data is also shown on the 
plot, with the number of standard deviations each observation devi- 
ates from the mean of the appropriate set of simulations. The distri- 
bution of this skewness statistic is not significantly altered between 
simulations that are consistent with WMAPl or WMAP3 data. The 
observed statistics for the WMAPl and WMAP3 data are similar 
but the slightly lower value for WMAP3 is now more apparent. 

To quantify the statistical significance of the detected devia- 
tion from Gaussianity we consider two techniques. The first tech- 
nique involves comparing the deviation of the observed statistic to 
all statistics computed from the simulations. This is a very conser- 
vative means of constructing significance levels. The second tech- 
nique involves performing a x' test. In both of these tests we re- 
late the observation to all test statistics originally computed, i.e. 
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Figure 2. Histograms of real Morlet wavelet coefficient skewness 
(an = 550'; y = 72°) obtained from 1000 Monte Carlo simulations. His- 
tograms are plotted for simulations in accordance with WMAPl (green) 
and WMAP3 (blue) observations. The observed statistics for the WMAPl 
and WMAP3 maps are shown by the green and blue lines respectively. The 
number of standard deviations these observations deviate from the mean of 
the appropriate set of simulations is also displayed. 




Figure 3. Histograms of normahsed x" test statistics computed from real 
Morlet wavelet coefficient statistics obtained from 1000 Monte Carlo simu- 
lations. Histograms are plotted for simulations in accordance with WMAPl 
(green) and WMAP3 (blue) observations. The x^ value computed for the 
WMAPl and WMAP3 maps are shown by the green and blue lines respec- 
tively. The significance of these observations, computed from the appropri- 
ate set of simulations, is also displayed. 



to both skewness and kurtosis statistics-. For a more thorough de- 
scription of these techniques see McEwen et al. (2005a). Searching 
through the 1000 WMAP3 simulations, 51 maps have an equiva- 
lent or greater deviation that the WMAP3 data in any single test 
statistic computed for that map. Using the very conservative first 
technique, the significance of the detection of non-Gaussianity in 
the WMAP3 data may therefore be quoted at 94.9 ± 0.7% (an ex- 
pression for the Icr errors quoted on significance levels is derived 
in Appendix A). The distribution oi x^ values obtained from the 
simulations is shown in Figure 3. The x^ value obtained for the 
data is also shown on the plot. The distribution of the x^ values 
is not significantly altered between simulations that are consistent 
with WMAPl or WMAP3 data. The x^ value computed for the 
data, however, is significantly lower for the WMAP3 data. Com- 
puting the significance of the detection of non-Gaussianity directly 
from the x^ distribution and observation, the significance of the de- 
tection of non-Gaussianity in the WMAP3 data may be quoted at 
97.2±0.5%. Using both of the techniques outlined above the signif- 
icance of the detection of non-Gaussianity made with the WMAP3 
data is slightly lower than that made with the WMAPl data. Nev- 
ertheless, the non-Gaussian signal is still present at a significant 
level. 

A wavelet analysis allows the spatial localisation of interest- 
ing signal characteristics. The most pronounced deviations from 
Gaussianity in the WMAP data may therefore be localised on the 
sky. The real Morlet wavelet coefficients of the WMAP3 data cor- 
responding to the most significant detection of non-Gaussianity on 
scale flu = 550' and orientation y = 72° are displayed in Fig- 
ure 4. Thresholded wavelet coefficient maps for both the WMAPl 
and WMAP3 data are also shown in order to localise the most 
pronounced deviations. The regions localised in the WMAPl and 
WMAP3 data are very similar, although the localised regions ap- 
pear slightly more pronounced, in the sense that the peaks are 
larger, in the WMAP3 data. To investigate the impact of localised 
regions on the initial detection of non-Gaussianity, the analysis is 

^ Although we recognise the distinction between skewness and kurtosis, 
there is no reason to partition the set of test statistics into skewness and 
kurtosis subsets. The full set of test statistics must be considered. 



repeated with the WMAP3 localised regions excluded from the 
analysis. The resulting skewness statistics are shown by the dashed 
line in Figure 1. Interestingly, the highly significant detections of 
non-Gaussianity are eliminated when these localised regions are 
removed. 

In our previous non-Gaussianity analysis (McEwen et al. 
2005a) we also performed a preliminary noise analysis and found 
that noise was not atypical in the localised regions that we detect. 
The localised regions have not changed in the WMAP3 data, hence 
we do not expect this finding to change. In an additional work of 
ours (McEwen et al. 2006) that investigated a Bianchi Vllh compo- 
nent as a possible source of non-Gaussianity - which, incidentally, 
we found not to be the predominant source of non-Gaussianity - we 
performed a preliminary analysis of foregrounds and systematics. 
We concluded that foregrounds or systematics were not the likely 
source of the detected non-Gaussianity. Again, we do not believe 
this finding to change in the WMAP3 data since both foregrounds 
and systematics are treated more thoroughly. 



4 SUMMARY AND DISCUSSION 

We have repeated on the WMAP3 data the directional spherical 
real Morlet wavelet analysis used to make a significant detection 
of non-Gaussianity in the WMAPl data (McEwen et al. 2005a). 
The non-Gaussian signal previously detected is indeed present in 
the WMAP3 data, although the significance of the detection is re- 
duced. Using our first very conservative method for constructing 
significance measures we find the significance of the detection of 
non-Gaussianity drops from 98.3 ± 0.4% to 94.9 ± 0.7%. Using our 
second technique for constructing significance measures, which is 
based on ax' analysis, the significance of the detection drops from 
99.3 ± 0.3% to 97.2 ± 0.5%. We have no intuitive explanation for 
this drop in significance. 

The most likely sources of non-Gaussianity were also lo- 
calised on the sky. We detect the same regions in the WMAP3 
data as found in the WMAPl data, although the localised regions 
extracted appear slightly more pronounced in the WMAP3 data. 
When all localised regions are excluded from the analysis the data 
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(a) WMAP3 wavelet coefficients 




(b) WMAP3 thresholded wavelet coefficients 




(c) WMAPl thresholded wavelet coefficients 

Figure 4. Real Morlet spherical wavelet coefficient maps and thresholded 
versions (an = 550'; y = 72°). To localise most likely deviations from 
Gaussianity on ffie sky, the coefficient map is thresholded so that only those 
coefficients above 3<t (in absolute value) remain. All sky maps (here and 
subsequently) are illustrated in Galactic coordinates, with the Galactic cen- 
ter in the middle. 



is consistent with Gaussianity. An interesting structure is extracted 
in the upper-left region of the thresholded maps (see Figure 5), i.e. 
in the vicinity of Galactic coordinates (Z, b) = (120°, 25°). In a fu- 
ture work we intend to use optimal filters on the sphere, in conjunc- 
tion with our fast CSWT analysis tool, to search for cosmic strings 
in the CMB, a possible source of the non-Gaussianity that we have 
detected in both the WMAPl and WMAP3 data. 
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Pr(x|p) = "C,p^(l-p)"--. 

The maximum likelihood (ML) estimate p of p is most easily given 
by maximiming the log-likelihood: 

51nPr(x|p) X 

= => p = - , 



dp 

which recovers the intuitive result. Approximating the shape of the 
likelihood near its peak by a Gaussian, we may approximate the 
standard deviation of p by 

1/2 



d-\nPr{x\p) 



dp'- 



x(n — x) 



APPENDIX A: ERRORS ON SIGNIFICANCE LEVELS 

In this appendix we derive the standard deviation of a significance 
level determined from Monte Carlo (MC) simulations. Suppose we 
perform n independent MC simulations. Let p denote the probabil- 
ity that an MC simulation chosen at random has a value for some 
test statistic that is larger that the corresponding value derived from 
the real data (hence p is the underlying significance we attempt 
to estimate). Choosing an MC simulation at random and defining 
whether it has a test statistic greater than that of the data thus cor- 
responds to a Bernoulli trial with a probability of success equal to 

P- 

Suppose we observe x successes in the n MC simulations. The 
likelihood for x is 
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